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Conditions required for the streaming effect and the optical-phonon transit-time resonance to 
take place in a compensated bulk GaN are analyzed in detail. Monte Carlo calculations of 
the high-frequency differential electron mobility are carried out. It is shown that the negative 
dynamic differential mobility can be realized in the terahertz frequency range, at low lattice 
temperatures of 30-77 K, and applied electric fields of 3-10 kV/cm. New manifestations of 
the streaming effect are revealed, namely, the anisotropy of the dynamic differential mobility 
and a specific behavior of the diffusion coefficient in the direction perpendicular to the applied 
electric field. The theory of terahertz radiation transmission through the structure with an 
epitaxial CaN layer is developed. Conditions for the amplification of electromagnetic waves in 
the frequency range of 0.5-2 THz are obtained. The polarization dependence of the radiation 
transmission coefficient through the structure in electric fields above 1 kV/cm is found. 
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1. Introduction 



The study of the streaming effect in semiconduc- 
tors was started in where it was demonstrated 
that a special streaming regime of electron trans- 
port can be realized in polar semiconductors at low 
temperatures and in strong enough applied electric 
fields. This regime of electron transport is charac- 
terized by the emergence of a specific quasiballistic 
motion of electrons in the momentum space. In other 
words, under the influence of a strong electric field, 
electrons are accelerated to the energy of an opti- 
cal phonon, hwo, almost without collisions. Having 
reached the energy ^ ^wq, the electron stops al- 
most instantly and emits a polar optical phonon. Af- 
terwards, a new cycle of acceleration begins. Such a 
cyclic motion of electrons results in the appearance 
of a strongly anisotropic streaming-like distribution 
function of electrons in the momentum space. The 
distribution function becomes strongly elongate along 
the electric field direction, being almost completely 
confined in the passive energy range, tp < fvjjQ. 
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For tlie electron streaming to be realized in semi- 
conductor materials, a number of conditions should 
be satisfied. 

(I) There must be < E On the 
one hand, the magnitudes of applied electric fields, 
_E, must reach the values close to the characteris- 
tic streaming field = po/e{Tp), where po is 
the electron momentum corresponding to the opti- 
cal phonon energy, and (tj,) is the averaged time of 
the electron momentum relaxation in the passive en- 
ergy range. On the other hand, the applied fields 
must be much lower than a certain characteristic 
field of the electron runaway effect, = Po/eT+, 
where is the emission time of a polar optical 
phonon. 

(II) For condition I to be satisfied, the mobility of 
electrons in low fields must be high enough, and, si- 
multaneously, the interaction between electrons and 
optical phonons must be strong, (rp) ^ Tq . 

(III) The lattice temperature T must be low, namely, 
fcbT < hujo^ where kh is the Boltzmann constant. 

(IV) The electron concentration must be low to 
avoid the electron-electron (e — e) scattering. At 
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least, the time of electron-electron collisions should 
exceed (r^), Te-e > (rp). 

The theoretical calculations of electric character- 
istics for hot electrons under streaming conditions 
meet difficulties in searching for a strongly nonequi- 
librium distribution function for charge carriers. The 
standard routine applied to the solution of the Boltz- 
mann kinetic equation, which uses the expansion of 
the distribution function in a series of spherical har- 
monics [2|,[3|, is not effective, being too cumbersome, 
because the harmonics of high orders (higher than the 
second one) have to be taken into account. For the 
same reason, the conventional diffusion approxima- 
tion [2| is unsuitable to describe the streaming. The 
electron temperature approximation j^] can be used 
only in the case of high electron concentrations, when 
the e — e interaction is the dominating mechanism of 
scattering. 

However, the approach proposed in [5] turned out 
to be rather fruitful and more adequate for the search 
of a strongly anisotropic distribution function. The 
Baraff method uses the approximation of a distri- 
bution function by the sum of isotropic and needle- 
shaped components. This method was widely applied 
in analytical researches of many types of problems 
dealing with the streaming effect fo'-Jq . Note that the 
Baraff method, in essence, gives rise to an approxi- 
mate solution of the Boltzmann kinetic equation and 
provides the exact solution only in the limit of perfect 
streaming, when E >• and = 0. 

Nowadays, the numerical methods got a wide popu- 
larity when being applied to the solution of the Boltz- 
mann kinetic equation. The most effective of them is 
the Monte Carlo method. It enables the exact solu- 
tion to be obtained for the Boltzmann kinetic equa- 
tion in a wide range of electric fields [E,J^]. With 
the help of the Monte Carlo method, it was shown 
that the streaming-like distribution function of elec- 
trons can be formed in polar semiconductor s in a dc 

^(str) 



11H13|. 



electric field with the amplitude E > 

Experimental confirmations of the formation of 
quasiballistic charge carrier motion in strong elec- 
tric fields were obtained in the 1980s, when study- 
ing the current-voltage characteristics of submicronic 
diodes fabricated on the basis of pure GaAs, InAs, 
and InP |l4ll5l |. In those experiments, the oscillatory 
dependence of the static differential conductivity on 
the applied electric bias was observed at low enough 
(helium) temperatures, with the period of oscillations 
corresponding to the value of Itwq / e. For A^^^B^ com- 
pounds, the oscillations of the static differential con- 



ductivity were observed in strong magnetic fields as 
well [l6|. 

Intensive researches of the streaming effect exe- 
cuted within the last decade involved a new class 
of wide-bandgap semiconductor materials, namely, 
group-Ill nitrides. Unique properties of nitride com- 
pounds (TtI such as, in particular, a large energy of 
polar optical phonon, a large value of Frohlich con- 
stant, and a relatively low effective mass (for GaN, 
those parameters are huo/k-b = 1000 K, ap = 0.4, and 
m* /me = 0.2, respectively) considerably improve the 
conditions for the streaming to take place. The Monte 
Carlo calculations carried out for GaN, InN, and AIN 
compounds f23-[26t showed that the streaming-like 
distribution function of electrons emerges at tempera- 
tures of 10-150 K and in fields of 1-30 kV/cm. It was 
demonstrated that the drift velocity of electrons, Vd, 
and their average energy, (e), saturate in this interval 
of fields and approach the values Vb/2 (Vb = pa/m*) 
and ^iwq/S, respectively. 

Modern researches of the streaming effect in ni- 
trides have a pronounced application aspect, namely, 
they are closely connected with the problem of de- 
veloping the terahertz radiation sources. It was 
shown theoretically (l8l . [l9| that the streaming regime 
can be accompanied by the emergence of a dynamic 
electric instability. This hypothesis was later con- 
firmed experimentally for InP at helium tempera- 
tures [23. The frequency dependence of the dy- 
namic (high-frequency) mobility in a system of 
streaming electrons has an oscillatory alternating- 
sign behavior. There exist the frequency intervals, 
in which Re[/itj] < 0. These frequency intervals are 
located near the characteristic transit-time frequency 
vr = eE/pQ. The latter corresponds to the reciprocal 
electron acceleration time in the static field E until 
the polar optical phonon energy is reached. The ap- 
pearance of the negative dynamic differential mobility 
is associated with the effect of electron bunching in 
the momentum space (2ll - |23t . 

Driven by an ac electric field with the resonance 
frequency w ^ttvh, the majority of electrons move 
in antiphase with the oscillations of this field, which 
results in the field strengthening. This effect is called 
the optical phonon transit-time resonance (OPTTR) 
effect or, shortly, the transit-time resonance. Its at- 
tractive feature consists in that the frequency and 
the amplitude of this resonance can be regulated by 
varying the strength of a dc electric field. This cir- 
cumstance opens wide perspectives for the creation 
of high-frequency sources of new types. 
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Calculations of the dynamic differential mobility 
in doped bulk nitrides fi?', '28'] and high-quality ni- 
tride quantum wells [,29. . .30i] showed that the dynamic 
negative differential mobility (DNDM) can reach sev- 
eral hundreds of cm^/V/s in the frequency range 
from 0.5 to a few terahertz, in electric fields of 1- 
10 kV/cm, and within the temperature interval of 
10-77 K. Similar conditions for the DNDM in GaN 



quantum wells were obtained in works j3lll32| . Note, 
however, that the e — e scattering was not taken into 
account in those calculations. The electron-electron 
coupling can change the conditions of existence for 
the DNDM very much. For a GaN quantum well, it 
was shown [ssl that the DNDM amplitude substan- 
tially decreases already at the electron concentration 
Tie = 10^^ cm~^ (the extrapolation of this value to 
the bulk sample gives the critical value of electron 



concentration n„ 



10 



16 



10" cm"3). At high elec- 



tron concentrations, when the e — e scattering dom- 
inates, the DNDM does not appear [sij. The nega- 
tive influence of the electron-electron scattering can 
be avoided with the help of the compensation of free 
carriers. A high compensation degree allows the elec- 
tron concentration to be reduced, hence making the 
electron-electron scattering not substantial. The case 
of compensated GaN was not discussed in detail in 
the literature. However, this case is important, be- 
cause it provides better conditions for the streaming 
and transit-time resonance effects to be observed. 

The main purpose of our researches was to reveal 
additional features in the electric characteristics of 
compensated GaN, which could definitely testify to 
the streaming emergence and could be identified in 
future experiments. In particular, we calculated, for 
the first time, the field dependence of the transverse 
diffusion coefficient and the frequency dependences of 
nonzero components of the dynamic mobility tensor 
in strong enough electric fields. It is of importance 
that those dependences can be observed in electro- 
gradient experiments and optical experiments deal- 
ing with the transmission of electromagnetic radia- 
tion with a given polarization through a structure 
that contains a layer of compensated GaN. For the 
calculations of stationary and high-frequency param- 
eters of the electron gas in strong electric fields, the 
Monte Carlo numerical method was applied. 

The paper is organized as follows. In Section [21 
the model of electron transport is described. In Sec- 
tion [Sj the features in the electron distribution func- 
tion that emerges in a constant electric field are dis- 
cussed, and the dependences of electric parameters 



of the electron gas on the amplitude of an applied 
field are analyzed. In Section IH the effect of transit- 
time resonance is studied, the spectra of the high- 
frequency mobility are presented for various relative 
orientations of dc and high-frequency electric fields, 
and the existence conditions for the DNDM are ana- 
lyzed. In Section [SJ the theory of terahertz radiation 
transmission through a structure with a thin epitax- 
ial GaN layer is developed. The main conclusions are 
summarized in Section |6l 

2. Model of Electron Transport 

Bulk GaN with a cubic modification and with the 
concentration of ionized impurities Ni = 10^^ cm~^ 
is considered. The concentration of electrons is sup- 
posed to be rie < Ni, i.e. the semiconductor is sug- 
gested to be compensated. The electron transport 
was simulated using the single-particle Monte Carlo 
method. The basic algorithms applied at the sim- 
ulation were standard; they are described in detail 
in works [l3|. In our model, we considered the 
processes of electron scattering by acoustic and po- 
lar optical phonons, as well as by ionized impurities. 
The dispersion law for electrons was assumed to be 
parabolic, and all the processes of electron scattering 
were supposed to occur only near the bottom of the 
lowest F- valley. Explicit expressions for the probabil- 
ities of electron scattering by acoustic and polar opti- 
cal phonons can be found in works |lol . l35| . The elec- 
tron scattering by ionized impurities was considered 
within the approach described in work (lo| : it was 
found to be more correct for compensated semicon- 
ductors in comparison with the conventional Brooks- 
Herring and Conwell-Weisskopf models. This ap- 
proach is described in Section r2 2. II in more details. 

2.1. Electron scattering by ionized impurities 

As a rule, the electron scattering by ionized impurities 
is considered in the framework of either the Brooks- 
Herring (BH) or the Conwell-Weisskopf (CW) model. 
In these models, the screening of the field, which is 
induced by impurity ions, by conduction electrons is 
taken into account in different ways. In particular, 
the BH model uses the screened Coulomb potential 



V{r) = — exp(-r/AD), 
Kor 



(1) 



where ZiC is the charge of an impurity ion, kq is the 
dielectric permittivity. Ad — {nok\,T / ATTe^ne)^^"^ is 
the Debye screening length, and rig is the electron 
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concentration. In the CW model, the unscreened 
Coulomb potential is cut off at the distance b = 
(3/47rA^j)^/'^ between impurity ions, and the mini- 
mum scattering angle ^min for an electron with energy 
e is determined by the formula 



ctg 



2ebKQ 



(2) 



Which of those models should be used depends on 
the ratio between Ad and b. In a semiconductor with 
a high compensation degree (rie <^ Ni), the inequal- 
ity Ad ^ & is obeyed so that the CW model proves 
to be more adequate for applications. For a heavily 
doped semiconductor, in which all impurities are ion- 
ized {Ni = He), the inverse inequality Ad "C & can 
be valid. In this case, the BH model is reasonable to 
be used. 

In the cases where Ad ^ b, it was suggested to 
use the CW model with the screened Coulomb poten- 
tial, rather than the "pure" one. In this model, the 
probability for an electron to transit from the initial 
state described by the wave vector k into the state 
with the wave vector k' within a unit time interval, 
provided that the electron is scattered by ionized im- 
purities, is given by the expression 



k,k 



Ad' + (k - k)- 



X(5(ek' -eic), 



(3) 



where is the normalization volume. In order to cal- 
culate the probability of the electron transition from 
the initial state into any other one within a unit time 
interval (the scattering rate), the quantity W^k.k' has 
to be multiplied by r2/(27r)'^, and the product has to 
be integrated over all k'-values, bearing in mind that 
the angle 9 between the vectors k and k' changes from 
^min to TT. For the parabolic dispersion law, we ob- 
tain the following formula describing the scattering 
rate for an electron with energy e: 



K2m*l/2fl/2 



1 



1 



Ed + 4e sin 



(4) 

where ed ~ ?i^/2to*Aq. Formula ^ gives the limit- 
ing transitions to the BH model as Oyaiu and to 
the CW one as ed — >■ 0. 

In a compensated GaN with the concentration of 
ionized impurities Ni = 10^^ cm~^, the concentra- 
tion of electrons rie — 10^^ cm"'^, and at a temper- 
ature of 30 K - below, it will be demonstrated that 



these parameters are the best for the realization of 
the OPTTR - the values 6 = 28 nm and Ad = 35 nm 
turn out close to each other so that it is formula (|4]) 
that should be used. 

2.2. Total scattering rate 

In our transport model, the total scattering rate is 
equal to the sum rtot = ''ac + nmp + ^op, where rimp, 
Tac, and Top are the probabilities of the electron scat- 
tering by ionized impurities, acoustic phonons, and 
polar optical phonons, respectively. The probabili- 



ties rac = r+ + Tac and r^n = r+ 



r^p take into ac- 



op ^ ■ op 

count the emission and absorption (— ) processes 
of acoustic and polar optical phonons, respectively. 
Note that, within the actual temperature interval, the 
electron scattering by acoustic phonons is almost elas- 
tic so that r+ ~ r~. At the same time, the mecha- 
nism of electron scattering by polar optical phonons is 
essentially inelastic. Therefore, the relation between 
the probabilities and r~p strongly depends on the 
lattice temperature and the electron energy. 

Figure 1 demonstrates the dependences of the total 
probability of electron scattering, rtot, on the electron 
energy in bulk GaN with the concentration of ionized 
impurities Ni — 10^^ cm^'^ and the concentration of 
electrons rig = 10^^ cm^'^ calculated for two lattice 
temperatures, T = 30 and 300 K. One can easily see 
a large difference between rtot-values in the passive 
(e < huJo) and active (e > huJo) energy regions. It 
can be explained by the fact that the electrons with 
the energy e < hujo are mostly scattered by ionized 
impurities and acoustic phonons, whereas the main 
scattering process for electrons with e > hujQ is driven 
by the more intensive spontaneous emission of polar 
optical phonons. For instance, for thermal electrons, 
i.e. electrons with the energy e = 0.028 x Hluq, the 
probabilities of their scattering by ionized impurities 



and acoustic phonons equal ' 



imp 



= 1.6 X 10^2 g-i and 



Tac = 1.6 X 10^^ s^^, respectively, at a lattice temper- 
ature of 30 K (curve 1 in Fig. 1). For an electron 
with the energy e = 1.2 x Hluq, the probability of the 
polar optical phonon emission is r+ = 4.5 x lO^'^ s^^ 



whereas r. 



10 



11 



and r\. 



7.5 X 10^2 s-\ At 



the same time, the polar optical phonon absorption 
is practically absent, r^^ ^ 10^ s^^. 

At room temperature of the lattice, T = 300 K 
(curve 2 in Fig. 1), the value of rtot for electrons with 
the energy e < hujQ turns out several times larger 
owing mainly to the growth of the role of inelastic 
processes in the optical phonon absorption. For ex- 
ample, for electrons with the energy e = 0.28 x huJo, 
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Fig. 1. Dependences of the total electron scattering rate 
on the electron energy at the temperature T = 30 (-Z) and 
300 K {2) 



we have rimp = 5.5 x 10^^ s ^, Tac = 5.1 x 10^^ s ^, 



and = 3.3 x 10 



" . At the same time, for 
an electron with the energy in the active range, e = 
1.2 X fvjjf)^ we have rac — 10^^ 



imp 



- 10^3 s-i, 



' op — " ^ , and 

The dependence of rtot on the electron energy, 
which was described above, is inherent to all nitrides. 
Note that the larger the difference between the rtot- 
values in the active and passive energy ranges, the 
more favorable are the conditions for the emergence 
of the streaming effect. 

3. Steady-State Electron Characteristics 

In the streaming regime, the electron motion becomes 
quasiperiodic, which finds its reflection in the distri- 
bution function of electrons and in the transport char- 
acteristics of the electron gas. The qualitative estima- 
tion of the characteristic electric field E'^^*''^^ , in which 
the streaming regime is realized, can be done on the 
basis of low-field mobility values for electrons. There- 
fore, instead of immediately proceeding to the analy- 
sis of stationary electron parameters, it is reasonable 
to discuss the dependence of the low-field electron 
mobility on the lattice temperature in detail. 

3.1. Low-field electron mobility 

At low temperatures (fcbT < hujo), when only the 
elastic mechanisms of scattering are actual, the elec- 
tron mobility can be calculated analytically with the 
use of the r-approximation. namely. 



MO = —{Tp), 
m* 

44 



(5) 



3V^(fcBr)5/: 



oo 



(6) 



is the statistically averaged momentum relaxation 
time. Here, is the inverse time of momentum re- 
laxation. According to Matthiessen's rule. 



(7) 



where Tp.ac and Tp^imp are the inverse momentum re- 
laxation time at acoustic phonons and ionized impu- 
rities, respectively. The explicit expressions for Tp^ac 
and rp^imp can be found in many manuals (see, e.g., 
textbook [36]). If the process of scattering by ion- 
ized impurities is examined in the framework of the 
Conwell-Weisskopf model with screened Coulomb po- 
tential, the following formula has to be used for the 
calculation of r„ 



' p.inip ■ 



nZfe^N, 1 



' p,imp 



4eeD(l - sin^ 



3/2 



log 



CD +4e 



CD + 4e sin 



2 e„ 







(eo +4e)(eD +4esin2 ^) 



(8) 



The dependence of Tp^mp on the electron energy 
e is shown in the inset in Fig. 2. At high energies, 
the value of rp^imp decreases with the growth of e as 
e~^/^. Such a behavior of rp. imp is explained by the 
fact that faster electrons are mainly scattered at small 
angles. For the sake of comparison, the probability 
of the electron scattering by ionized impurities rimp 
grows with the electron energy (see Fig. 1). 

Note that, at high temperatures, when the role of 
the inelastic scattering mechanisms becomes essen- 
tial, the T-approximation loses its meaning, so that 
expression ([5]) cannot be used. For the exact calcu- 
lation of the electron mobility in a low field, /io, in a 
wide temperature range, either of the following proce- 
dures can be applied: (1) the Monte Carlo method is 
used to calculate the dependence of electron drift ve- 
locity Vdr in the field with strength E: afterwards, the 
electron mobility can be determined from the slope of 
the curve Vdr(£') [111; or (2) the Monte Carlo method 
is used to calculate the diffusion coefficient Dq; after- 
wards, the electron mobility can be determined from 
the Einstein relation /io = eDo/khT [38]. In the case 
of low fields, the second way turns out to be more 
accurate and less dependent on the statistical noise 
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produced by the Monte Carlo calculations. There- 
fore, it was the second way that was selected by us to 
calculate the dependence of the electron mobility in 
low fields on the lattice temperature. 

As is seen from Fig. 2, the dependence of the low- 
field mobility /zo on T is nonmonotonous. In the tem- 
perature interval 30-120 K, the /io- values increase 
with T, because electrons are mainly scattered by 
acoustic phonons and ionized impurities, with the 
latter process dominating. It is of interest that the 
growth of fiQ at higher T is associated with the de- 
creasing dependence of Vp, 

imp on the energy. As the 
lattice temperature increases, the fraction of high- 
energy electrons, for which scattering by ionized im- 
purities is less intensive, grows. In Fig. 2, the results 
of calculations of the low-field mobility carried out ei- 
ther in the framework of the Monte Carlo method or 
with the use of formulas ([S])-® are evidently identi- 
cal within the temperature interval T = 30 — 120 K. 
Starting from a temperature of 120 K, the electron 
mobility diminishes with the growth of T, which 
is connected with the enhancement of the role of the 
mechanisms of scattering by acoustic and polar opti- 
cal phonons. If the temperature continues to grow, 
the role of the electron scattering by polar optical 
phonons increases considerably. In Fig. 2, this fact is 
illustrated by an increasing discrepancy between the 
exact mobility value obtained by the Monte Carlo 
method and its approximation calculated by formu- 
las (O-®. At room temperature, the mobility of 
electrons is mainly governed by their scattering by 
polar optical phonons. The values calculated by us 
are close to those measured in high-quality GaN epi- 
taxial layers grown up on AI2O3 substrates [s^. 

At T = 30 K, the electron mobility equals 
5000 cm^/V/s so that the characteristic electric field 
of the streaming is E^^'^ = 8 kV/cm; whereas the elec- 
tron mobility at T = 77 K is close to 10000 cm^/V/s, 
and E""^' = 4 kV/cm. 

3.2. Electron distribution function 

Generally speaking, the distribution function of elec- 
trons in a bulk semiconductor is a function of three 
momentum variables, F[px,Py,Pz)- However, the dis- 
tribution function of electrons in a uniform dc electric 
field has the axial symmetry with respect to the field 
direction. Hence, if an electric field directed along the 
z-axis is applied to the semiconductor, it is enough to 
analyze the distribution functions of electrons in the 
momentum space in two directions, namely, along the 



> 

S 




4000 



T, K 

Fig. 2. Dependences of the low-field electron mobility on 
the lattice temperature calculated exactly by the Monte Carlo 
method (solid curve) and approximately by formulas ll5t— lISll 
(dashed curve). Ni = 10^^ cm""^ and rie = 10^^ cm~'^ 



field, f{pz), and across it, f{px)- The distribution 
functions f{pz) and f{px) are introduced as follows: 



f{Pz) = I I dpxdpyF{px,py,pz)/N 
and 

f{Px) =11 dpzdpyF{px,py,pz)/N, 
where 



N = 



dpxdpydpzF{px,Py,Pz) 



is the normalization integral. 

The calculated distributions of electrons in the mo- 
mentum space along and across the applied electric 
field are shown in Fig. 3 (panels a and &, respectively). 
Already in a field of 0.5 kV/cm, the system of elec- 
trons becomes nonequilibrium, and the difference be- 
tween the electron distributions along and across the 
electric field is well distinguished. The distribution 
of electrons along the electric field has an apprecia- 
ble asymmetry, which is associated with the existence 
of two groups of electrons, namely, low-energy elec- 
trons, which are isotropically distributed in the space 
of their momenta, and high-energy ones with the mo- 
menta directed along the electric field. As the electric 
field grows, the number of such monodirected elec- 
trons increases, and the electron distribution becomes 
more and more anisotropic. In the perfect streaming 
limit, all electrons have their momenta directed along 
the field, and the distribution function of electrons 
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P/Po 

Fig. 3. Distributions of electrons in the momentum space 
(a) along and (6) across the electric field in GaN. Ni = 
= IQis cm-3, ne = IQis cm-3, and T = 30 K 




u.u^ , , , 1 , 1 . 1 , 1 

2 4 6 8 10 

E, kV/cm 

Fig. 4- Dependences of (a) drift velocity and (b) average 
energy of electrons on the electric field amplitude. Vb = 
4 X 10 cm/s is the characteristic velocity. The other parame- 
ters of the material are the same as in Fig. 3. Points correspond 
to the field values, at which the distribution functions depicted 
in Fig. 3 were calculated 

along the field has a step-like shape. From Fig. 3, a, 
it is evident that, for the fields of 3-8 kV/cm, the 
shape of the electron distribution function along the 
field is close to the step-like one. At the same time, 
the distribution function of electrons across the field, 
fipx), remains symmetric. In the directions perpen- 
dicular to the field, the dominating fraction of elec- 
trons have momenta Px,Py < Po- For the growing 
field, one can observe a weak compression of tails of 
the transverse component of the distribution function 
f{px) (see Fig. 3,6). 



In the fields higher than 10 kV/cm, the streaming 
effect becomes destroyed, because electrons penetrate 
deeper into the active energy range. 

3.3. Steady-state electric characteristics 

The streaming regime can manifest itself through 
the emergence of characteristic features in the de- 
pendences of stationary electric parameters on the 
applied field. Among those, there are weak depen- 
dences of the drift velocity Va and the average en- 
ergy (e) on the applied electric field strength E. In 
Fig. 4, one can easily see that, in the range of low, 
i.e. pre-streaming, electric fields 0.1-1 kV/cm, the 
dependence of the drift velocity on the field is ap- 
proximately linear, but becomes sublincar if the field 
continues to grow. In the fields 3-10 kV/cm, the 
streaming type of electron transport is formed: the 
drift velocity practically does not change as the field 
grows, being close to Vb/2, which corresponds to the 
drift velocity in the perfect streaming model. A simi- 
lar scenario is observed for the field dependence of the 
average energy. In the developed streaming regime, 
the value of (e)(-E) approaches that of huJo/3, which 
is characteristic of the perfect streaming limit. It is 
worth noting that the dependences Vd{E) and {e){E) 
qualitatively do not differ strongly from each other 
within the temperature interval 30-150 K. 

At room temperature, the situation is different. 
The streaming regime does not emerge owing to the 
presence of the strongly inelastic scattering in the pas- 
sive energy range, the mechanism of which consists 
in the absorption of a polar optical phonon. In the 
fields 3-10 kV/cm, the dependence of the drift ve- 
locity on the field remains almost linear, the average 
energy almost does not change and remains close to 
the equilibrium value 3/2 x k^T. 

The emergence of the streaming regime can also 
be clearly traced by detecting a noiimonotonous field 
dependence of the transverse component of the aver- 
age electron energy, (e^) = ((p-J. +Py)/2rn*), where 
P:c and py are the electron momenta in the directions 
perpendicular to the field. In Fig. 5, the field depen- 
dences of {e±) calculated for various lattice temper- 
atures are depicted. At cryogenic temperatures and 
in low fields (less than 1 kV/cm), the value of {e±) 
grows owing to the heating of the electron gas, which 
mainly manifests itself in the isotropic broadening of 
the electron distribution function in the momentum 
space. If the field amplitude increases further, {e±) 
decreases and saturates in fields that correspond to 
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the developed streaming, which is associated with the 
narrowing of high-energy tails of the electron distri- 
bution function in the directions perpendicular to the 
field (see Fig. 3,6). At room temperature, both (e^) 
and the total average energy (e) practically do not 
change as the field grows. 

Another attribute of the fact that the electron 
transport is in the streaming regime consists in a 
specific behavior of the field dependence of the elec- 
tron diffusion coefficient in the coordinate space. As 
was shown above, the electron distribution function 
across the field gets narrowed at the streaming for- 
mation. As a result, the diffusion motion of electrons 
across the field must be inhibited, which would lead 
to a reduction of the diffusion coefficient D± in the 
direction perpendicular to the field. This feature was 
discussed in work fs?!, when analyzing the stream- 
ing effect for the two-dimensional electron gas in the 
framework of the Baraff approximation. For today, 
a lot of researches have been carried out concerning 
the diffusion coefficient in bulk nitride samples [4^ in 
strong fields (up to 0.5 MV/cm). At the same time, 
the interval of moderate electric fields, at which the 
streaming effect becomes possible, was not given a 
proper attention. 

In Fig. 6, the field dependences D±{E) calculated 
for various lattice temperatures by the Monte Carlo 
method are shown. In heating fields of 0-0.25 kV/cm 
and at cryogenic lattice temperatures, the magnitude 
of D± drastically grows. If the field amplitude grows 
further and the streaming-like distribution of elec- 
trons is formed, the value of D± decreases. Starting 
from a field of 3 kV/cm, i.e. in the developed stream- 
ing regime, D± saturates and approached values of 
20-25 cm^/s, which are close to or even lower than 
the values for the equilibrium diffusion coefficients — 
the corresponding values are 13, 57, and 120 cm^/s 
at lattice temperatures of 30, 77, and 150 K. respec- 
tively. At room temperature, D± is almost indepen- 
dent of the field amplitude, varying from 50 cm^/s in 
the zero field to 30 cm^/s in a field of 10 kV/cm. 

Thus, a number of features in the field dependences 
of electric parameters, which point to the emergence 
of the streaming effect, are observed for the electron 
gas in a compensated bulk GaN with the concen- 
tration of ionized impurities Ni — 10^^ cm^'^, the 
electron concentration rig = 10^^ cm~^, in the fields 
of a few kV/cm, and at lattice temperatures of 30- 
150 K. It is worth noting that the streaming regime 
can be identified according to the results of exper- 
imental measurements of the drift velocity and the 
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Fig. 5. Transverse component of the average electron energy 
as a function of the field 



250 



200- 



E 
o 




E, kV/cm 

Fig. 6. Diffusion coefficient in the direction perpendicular to 
the field as a function of this field. The other parameters of 
the material are the same as in Fig. 3 



diffusion coefficient of electrons in the direction per- 
pendicular to the field. 

4. High-Frequency Electric Characteristics 

Surely, the most remarkable manifestation of the 
streaming phenomenon is the transit-time resonance 
effect, which is connected with the emergence of the 
negative dynamic differential conductivity as a re- 
sponse to an external high-frequency signal. In this 
section, we will discuss the spectra of the dynamic 
differential mobility (DDM), fj,^, calculated by the 
Monte Carlo method using the single-particle algo- 
rithm [1^,13, [4^. We will analyze comprehensively 
the conditions of existence for the negative dynamic 
differential mobility (DNDM) in a compensated GaN 
at various lattice temperatures and various relative 



ISSN 2071-0186. Ukr. J. Phys. 2013. Vol. 58, No. 1 



47 



G.I. Syngaivska, V.V. Korotyeyev 



TQOOO 




t -i — . — 1 — ■ — 1 — ■ — I I 

Q.& 0,2 0.4 0.6 0.8 1.0 
«/27i, THz 



Fig. 7. Dependence of DDM on the frequency calculated in 
the framework of the Drude-Lorentz model (thin curves) and 
by the Monte Carlo method (bold curves). Solid curves cor- 
respond to Re[/iuj] and dashed ones to Im[/ia;]. The dc field 
amplitude E = 0.5 kV/cm, T = 30 K, AT^ = 10^® cm-^, and 
He = lO^^ cm-3 




ao 0.5 i,p ij 



co/27c,THz 

Fig. 8. Dependences Re[/iaj] (solid curves) and Im[/ic^] (dashed 
curves) at T = 30 K and E = 1.5 (1) and 3 kV/cm (2). The 
spectra Re[/i;^] and Im[/i;^] at E = 3 kV/cm and T = 300 K 
axe shown in the inset 

orientations of dc and ac electric fields. The analysis 
of flu is started from the case where a dc field E and 
a low sinusoidal field with amplitude -C E, which 
are parallel to each other, are applied to the sample. 

4-1- Dynamic differential mobility 

In Fig. 7, the bold curves correspond to the frequency 
dependences of calculated for the low station- 
ary electric field E = 0.5 kV/cm using the Monte 
Carlo method. For comparison, the same depen- 
dences calculated in the framework of the Drude- 



Lorentz model, 

/-ic^ = Mo,b/(1 - iu!m*fio,E/e), 

where /Lto,B = dVd{E)/dE is the zero- frequency differ- 
ential mobility, are also shown by thin curves. In the 
field E = 0.5 kV/cm, the electron distribution func- 
tion still remains quasiisotropic, and, as is seen from 
Fig. 7, the shape of the /Zc^-spectrum is close to that 
of the Drude-Lorcntz one. 

If the dc field amplitude increases, the streaming 
regime starts to be formed, which manifests itself in 
the oscillatory dependence of the real and imaginary 
parts of DDM on the frequency. In Fig. 8, the /x^^- 
spectra calculated for E = 1.5 and 3 kV/cm are de- 
picted. Already at E ~ 1.5 kV/cm, the real part of 
DDM has a series of minima at the transit-time res- 
onance frequency and its higher harmonics, but still 
remains positive. At = 3 kV/cm. the real part 
of DDM changes its sign and becomes negative in a 
vicinity of the fundamental frequency of the transit- 
time resonance wtr/27r w 0.6 THz. Note that, at 
room temperature and in the same electric field, the 
real and imaginary parts of DDM do not reveal any 
features, being well described by the Drude-Lorentz 
model. In the inset of Fig. 8, bold (calculation by the 
Monte Carlo method) and thin (the Drude-Lorentz 
model) curves are close to each other. 

Hence, the amplitude and frequency windows of 
DNDM directly depend on the applied field magni- 
tude and the temperature of a sample. In the fre- 
quency range, where Re[/Lta;] < 0, the amplification of 
a high-frequency signal occurs, which is proportional 
to the DNDM amplitude. Therefore, it is expedient 
to determine the intervals for temperatures and dc 
electric fields, at which the DNDM takes place and 
reaches minimum negative values. 

4-2. Dependences of the OPTTR effect on the 
temperature and the field amplitude 

In Fig. 9, the frequency dependences of Re[/iij] are 
depicted only for the first actual resonance mini- 
mum. At higher-order minima, which are multi- 
ples of the transit-time frequency, the DNDM does 
not emerge. The most pronounced OPTTR effect 
is observed at a temperature of 30 K. The DNDM 
appears at E ^ 2 kV/cm and reaches a minimum 
value of —250 cm^/V/s in a vicinity of the frequency 
w/27r w 0.64 THz at w 3 kV/cm. As the dc field 
amplitude increases, the frequency windows of the 
DNDM shift toward the high-frequency range and are 



48 



ISSN 2071-0186. Ukr. J. Phys. 2013. Vol. 58, No. 1 



Electrical and High- Frequency Properties of Compensated GaN 



broaden out. At the same time, however, the DNDM 
amphtude decreases. In Fig. 9, it is clearly seen that, 
at T = 30 K, the DNDM appears in the field range 
within the limits 2 10 kV/cm and in the frequency 
interval from 0.38 to 2.5 THz (this interval is con- 
fined by a dash-dotted curve). However, if the field 
amplitude exceeds 10 kV/cm, a sufficient number of 
electrons can penetrate into the active range to vio- 
late the coherent motion of the majority of electrons, 
so that the DDM values become positive. 

At T = 77 K, the frequency windows with the 
DNDM still exist. However, the DNDM manifests it- 
self much weaker than at T = 30 K. The frequency in- 
terval, in which the DDM becomes negative, is much 
narrower at T = 77 K than at T = 30 K. At the 
fields within the interval 2.5-9 kV/cm, the DNDM 
emerges at frequencies between approximately 0.56 
and 2.2 THz. The largest DNDM amphtude is real- 
ized &t E = A kV/cm at a frequency of 0.9 THz and 
with the minimum Re[/i[j] « —100 cm^/V/s. 

At T = 150 K, the DDM also demonstrates the 
oscillatory dependence on the frequency; however, 
it does not become negative. It is of interest that, 
at this temperature, the low-field mobility /xq = 
9500 cm^/V/s (see Fig. 2), which is almost twice as 
large as that at T = 30 K (/io = 5000 cmVV/s). 
Hence, the conditions for the streaming effect to take 
place seem to be better at T = 150 K rather than 
at T = 30 K, and, consequently, the OPTTR effect 
should have manifested itself more strongly just in 
the former case. However, the opposite situation is 
actually observed. Such a disagreement can be ex- 
plained by the sensitivity of the OPTTR effect to the 
initial broadening of the equilibrium electron distri- 
bution function. For the distribution function with 
an anisotropy, which would be sufficient for the ap- 
pearance of the DNDM at T = 150 K, to be formed, a 
field stronger than that required at 30-77 K has to be 
applied to the semiconductor. However, in stronger 
fields, the DNDM does not appear, because the pen- 
etration of electrons into the active energy range be- 
comes substantial. 

4..3. Anisotropy of dynamic differential 
mobility 

In Section 4.2, the dependences of the real part of 
DDM on the frequency of the ac electric field and the 
amplitude of the dc one, which are parallel to each 
other, were discussed in detail. Under the stream- 
ing conditions, when the distribution function of elec- 
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Fig. 9. R,c[^tt^.]-spcctra in a vicinity of OPTTR frequencies 
at £ = 3 (1), 5 (2), 7 (,?), and 9 W/cm (4). The dash-dotted 
curve is the envelope of Re[jUu)]-niinima. Ni = 10^^ cm"** and 
ne = 10^^ cm-3 
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Fig. 10. Spectra of /iji (bold curves) and fi^ (thin curves) 
at T = 30 K and E = 1 (a) and 3 W/cin (6). The real 
and imaginary parts of DDM are shown by sohd and dashed, 
respectively, curves. The same spectra but at T = 300 K and 
E = 3 kV/cm are shown in the inset 

trons is anisotropic, one may expect that the electron 
response should depend on the relative orientation of 
the dc field and the varying signal. 

In Fig. 10, the dependences of DDM on the fre- 
quency of ac signals polarized in parallel and per- 
pendicularly to the direction of the dc electric field 
are exhibited. The quantities //i- and /i^fr describe 
the response of the electron system to the ac sig- 
nal with the parallel and perpendicular, respectively, 
polarizations. At low frequencies, the ^li- and fx^- 
values substantially differ from each other both in 
low fields (< 1 kV/cm), in which the distribution 
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Fig. 1 1 . Dependences /jjj, (bold curves) and fi-^ (thin curves) 
calculated from stationary characteristics Vd(E), (ej^ }(£?), and 
D±(E). iJq (squares) and /1q- (circles) values were calculated 
by the Monte Carlo method. The lattice temperature is 30 K 

function of electrons is characterized by a moderate 
anisotropy, and in the developed streaming regime 
= 3 kV/cm). 

The difference between ^i, and /i^^ survives even 
at the zero frequency, which is associated with the 
anisotropy of the electron distribution in the momen- 
tum space. In Fig. 11, the dependences of parallel, /ijl, 
and perpendicular, ^q, components of the differential 
mobility on the dc field strength E are shown. These 
dependences were obtained with the use of the Monte 
Carlo method, in which the /^jj- and /ij^-values were 
calculated for an extremely low frequency. There is 
another approximate method for the determination of 
those quantities, which uses the values of stationary 
parameters. In a given dc field E, « dVd{E)/dE. 
This relation is exact in the limit of low-signal re- 
sponse {Ei^ <^ E). To calculate /i^, the generalized 
Einstein relation can be used. 



with regard for the dependences {e±){E) (see Fig. 5) 
and D±{E) (see Fig. 6). The values calculated for 
/ijl and by the Monte Carlo method coincide with 
the corresponding values obtained from the station- 
ary characteristics Vd{E), {e±){E), D±{E), and for- 
mula ^ (see Fig. 11). The fact that both determina- 
tion methods bring about identical values for iiq{E) 
can be explained by an almost Maxwellian distribu- 
tion of electrons in the direction perpendicular to the 
dc field. 



From Fig. 11, one can easily see that the de- 
pendence ^J.Q{E) decreases more rapidly than the 
IJ.q{E) one with the growth of the dc field amplitude. 
Such a behavior of fJ-^iE) can be explained by the 
fact that the sublinearity manifests itself much more 
strongly in the dependence Vd{E), rather than in the 
D^{E)/{e±){E) one. In fields of 3-8 kV/cm, which 
correspond to the developed streaming regime, a con- 
siderable difference between the ^q{E)- and ^q{E)- 
values is retained. Hence, the experimentally ob- 
served substantial anisotropy of the differential mobil- 
ity in strong electric fields can serve as an additional 
proof of the streaming formation. 

In fields of 1-3 kV/cm, the distribution function of 
electrons in the momentum space remains symmetric 
in the direction perpendicular to the field. Conse- 
quently, the frequency dependence of /i^ is well de- 
scribed by the Drude-Lorentz model and does not re- 
veal an oscillatory behavior typical of /xli (see Fig. 10). 
It is worth to note that there is no substantial dif- 
ference between /i^- and /xJi-values in a dc field of 
1 kV/cm, except for at very low frequencies. Un- 
der the streaming conditions (e.g., at 3 kV/cm), the 
response of the electron system in the parallel con- 
figuration of the fields E and E^ differs cardinally 
from that in the perpendicular configuration. In the 
frequency interval, where becomes negative and 
/i^ remains positive, an ac signal with the polariza- 
tion perpendicular to the dc field, instead of being 
amplified, is effectively absorbed. This effect can be 
observed in optical experiments on the transmission 
of ultrahigh-frequency (terahertz) radiation through 
semiconductor structures. If a strong enough dc elec- 
tric field is applied to a sample, the latter is char- 
acterized by the anisotropy of DDM and acts as a 
polarizer for a non-polarized beam. The efficiency of 
such a polarizer depends on the magnitude of applied 
dc electric field and the sample temperature. For in- 
stance, aX T — 300 K, when the electron distribution 
function is still isotropic and the /i^- and /i|i-values 
practically coincide, the sample does not function as 
a polarizer. 

Recently, the collaborators of the the Terahertz 
laboratory at the Montpellier University carried out 
an experiment, in which they tried to register the 
OPTTR effect and find the DNDM by detecting 
the amplification of terahertz radiation transmitted 
through a heterostructure with GaN [4^. Unfortu- 
nately, reliable confirmations for the OPTTR effect 
were not obtained. To elucidate what one could ex- 
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pect of such experiments, we developed a theory de- 
scribing the hght transmission through a sample with 
an active epitaxial layer of compensated GaN in the 
OPTTR regime. 

5. Transmission of High-Frequency Radiation 
through a GaN Structure 

In modern experiments, the complicated multilayered 
structures grown up on dielectric substrates are used. 
As a rule, the thickness of the active zone, a thin GaN 
layer, has an order of a few micrometers, which is 
much less than the wavelength Ao of terahertz electro- 
magnetic radiation in vacuum. In similar structures, 
the thickness of dielectric substrate is much thicker 
than that of GaN layer, being, as a rule, of the same 
order of magnitude as Aq. In this section, we expound 
the theory of high-frequency radiation transmission 
through such structures. While studying the spec- 
tra of the transmission, reflection, and absorption co- 
efficients for high-frequency radiation, the frequency 
dependences of DDM were applied, which were cal- 
culated with the use of the Monte Carlo method (see 
Section g]). 

5. 1 . Theory of high-frequency radiation 
transmission through GaN 

Let high-frequency radiation pass through a struc- 
ture consisting of a dielectric substrate with the thick- 
ness ds and the dielectric constant k^, covered with 
a delta-like GaN layer. The GaN layer is charac- 
terized by the two-dimensional complex conductivity 
cr* = (j'^d + id{a'^ - Kow/47r), where a'^ = eneRe[/x^], 
(7^ = enelm[/x^], and d is the thickness of GaN 
layer. This expression takes the displacement cur- 
rent in the GaN layer into account. Let a plane wave 
Ep^i^(y) exp{—iujt) with the amplitude Ep^i^{y) and 
the frequency lu fall normally on the structure sur- 
face. The electric field Ep,i^{y) of this wave satisfies 
the Maxwell equations. 



d^E, 



P,UJ 



dy'^ 



75- 



V<0 
0<y<ds 

y>d. 



> E, 



'p,uj 



4:TTiujcr*;P 



Ep^ujS{y), 



(10) 



where the subscript p specifies the wave polarization 
either along {p — |j) or across {p — 1.) the dc electric 
field. In Eqs. ([TU)) . the whole structure is supposed to 
be in vacuum. The solution of system (|10p looks like 

Ep,u:{y) = 

ISSN 2071-0186. Ukr. J. Phys. 2013. Vol. 58, No. 1 



' Ap^^ eyiY>[ikoy) + Bp^^ exp(-ifeoJ/), y < 0, 
Cp^uj exp(ifcsy) -I- Dp^^ ex.p{-iksy), < y < dg, 
^Fp,^exp(ifcot/), y>ds, 

(11) 

where fcg = ^/c and kg = io^/Kljc are the wave num- 
bers of plane waves in vacuum and the substrate, re- 
spectively. The coefficients Ap^^, i?p,tj, Cp^^)-, ^p.tj, 
and Fp,^ are determined from the following condi- 
tions at the coordinate planes y = and y — dg'. 

Ep^^i-0)=Ep^U+0), 



^(-0)-^^(+0) = l^i?,.(0), 
dy dy c/ 



Ep^Uds-0)=Ep^^{dg+0), 



'^^P-^-{dg~0)^^^{ds-^0). 



dy 



dy 



(12) 



After standard transformations, we obtain the fol- 
lowing formulas for the transmission, Tp ^^ , and refiec- 



tion, Rp,uj, coefficients: 



T = 



1 



F' 



2 p//2 



2 7 4 

(/..-i)(^.-(i + r;j^-r-i) 



X sin^(fcsds) — "^^^^ ^ s[T[i(2ksds) 



4kc 



(13) 



R 



+ 



(^^.-i)(^.-(i-r;, j^-r;^) . 

— — sm [kgds 



4k., 



4k., 



sm{2ksds) 



(14) 



The absorption coefficient can be calculated by the 



formula Lp^^^ = 1 - Tp^^^ - Rp^^^ or 



T T^f I ^S 1 . 2 



sin^(fc^4) X Tp,, 



(15) 



In all Eqs. HI])- US]), the notations = 
= 47rRe [cr*'P]/c and r^ ,^ = 47rIm[(T*^P]/c are used. 
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Fig. 12. Spectra of the (a) transmission, (b) reflection, and 
(c) absorption coeflicients for the parallel (solid curves) and 
perpendicular (dashed curves) field configurations. Tguj,, iJsubi 
and -Lsub spectra are shown by dash-dotted curves. The param- 



eters of the GaN layer are Ni ■- 
d = 10-3 cm, Ko = 8.9, and E ■- 



10^" cm-3. Tie = cm--^, 
3 kV/cm. The substrate (sap- 
phire is supposed) parameters are ds = 0.04 cm and Ks = 12 



3 4 

E, kV/cm 

Fig. 14. Dependences of (a) the real (bold curves) and imag- 
inary (thin curves) parts of (solid curves) and ^j^ (dashed 
curves), (6) the transmission coefficients along (Ty bold 
curve) and across (Tx b, thin curve) the applied field, and 
(c) the polarization degree Pe on the field amplitude at a fre- 
quency of 0.64 THz. The parameters of the substrate and the 
active layer are the same as in Fig. 12 
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Fig. 13. (a) Transmission coefficient spectra for a single sub- 
strate and the (fc) absorption coefficient spectra for the struc- 
ture at ds = 0.036 (dash-dotted curves) and 0.04 cm (solid 
curves). The other parameters of the material are the same as 
in Fig. 12 

In the absence of a GaN layer, i.e. if there is only 
the substrate, the quantities ^ and T'^ ^ in formu- 
las ([T5| - (fT5)) equal zero, and the transmission, reflec- 
tion, and absorption coeflicients are determined in a 
standard way, as for a single dielectric wafer, namely. 



l + ^sm2(W 
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Lsnh = 0. (16) 

5.2. Spectra of the transmission, reflection, 
and absorption coefficients 

The transmission, reflection, and absorption coeffi- 
cients for high-frequency radiation (Fig. 12) were cal- 
culated in such frequency intervals and at such val- 
ues of dc electric field, when the DNDM manifests 
itself most strongly (Fig. 10,6). If an electromagnetic 
wave passes through a single substrate, then, for the 
given thickness and the dielectric constant of the lat- 
ter, we obtain a set of frequencies uj^ {r — 1, 2, 3, ...), 
which correspond to the Fabry-Perot modes of a 
plane-parallel dielectric wafer. Their values are given 
by the expression ui-rjlix — cr/{2dgyfKl) and corre- 
spond to Tsub = 1- The presence of a thin active 
GaN layer with a low concentration of electrons on 
the substrate slightly modifies the transmission coef- 
ficient spectra of the system. The parameters of the 
substrate and the active element were so chosen that 
the DNDM frequency window should coincide with 
any of the frequencies Ur- In the frequency interval 
0.6-0.7 THz, the wave with the polarization along the 
dc field, when passing through the sample, becomes 
amplified. In this frequency interval, the absorption 
coefficient becomes negative and reaches a minimum 
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value of -1.5% at about 0.64 THz (Fig. 12). The 
negative absorption coefficient means that the sum 
of intensities for passed and reflected waves exceeds 
the incident wave intensity. Hence, we may say about 
the amplification of the electromagnetic field by the 
active element. For the sake of comparison, the coeffi- 
cient of electromagnetic wave losses at this frequency 
under the wave reflection from a perfect silver mirror 
equals —0.5%. Therefore, the criterion of terahertz 
mode excitation in a resonator system composed of 
two plane-parallel metallic mirrors can be satisfied 
despite that the gain factor in the active layer is low. 

For a wave with the polarization perpendicular to 
the dc field, the absorption coefficient is positive, and, 
hence, there is no amplification of the electromag- 
netic field. 

It should be noted that the coefficient of electro- 
magnetic wave absorption in the subterahertz fre- 
quency range depends on the substrate parameters. 
In Fig. 13, the frequency dependences of L^^ are shown 
for two samples with identical parameters of their 
active elements, but diflFerent substrate thicknesses. 
If the substrate thickness dg = 0.036 cm, the fre- 
quency window of DNDM does not coincide with any 
of the frequencies Ur of Fabry-Perot modes. At such 
non- resonance parameters, the values of \L^\ and, 
hence, the amplification of the electromagnetic field 
are substantially lower than those in the resonance 
case where ds = 0.04 cm. Therefore, we may say 
about a selective role of the substrate in the amplifi- 
cation of the electromagnetic field in the subterahertz 
frequency range. 

In experiments dealing with the transmission of ra- 
diation through semiconductor structures, it is much 
more convenient to measure the transmission coeffi- 
cient at a given frequency by varying the amplitude 
of the applied dc field. In Fig. 14, a, the dependences 
of DDM on the applied field amplitude obtained in 
the parallel, and transverse, ^j^, field configura- 
tions at a frequency of 0.64 THz are shown. The real 
and imaginary parts of /J^ have an oscillatory behav- 
ior, and the DNDM is realized in a narrow interval 
of dc fields with the amplitudes of about 3 kV/cm. 
At the same time, the real and imaginary parts of 
almost do not change at that. At a frequency of 

0.64 THz, a substantial difference between and 
fi-^ is observed starting from the applied field values 
of 2-3 kV/cm. 

Experimentally, it is possible to observe the field- 
induced difference between the transmission coeffi- 



cients for electromagnetic waves with polarizations 
along, Tii.Bj Siiid across, T±_e, the dc field. In 
Fig. 14,6, the bold and thin curves illustrate the 
field dependence for T|| ^ and T±^e, respectively. A 
monochromatic beam, which initially was not polar- 
ized, after having passed through the sample, became 
partially polarized. The polarization degree of such 
a beam, Pe- depends on the dc field amplitude, as 
is shown in Fig. 14, c. The quantity Pe is defined as 
follows: Pe = iT\\,E-T±^E)/m^E + T^.E). From 
Fig. 14, c, one can see that the behavior of Pe repro- 
duces the oscillatory dependence of DDM ij}^ on the 
field, which points to the formation of the streaming 
and the appearance of the OPTTR. Beyond the range 
of resonance fields, the value of Pe monotonously 
grows, and the oscillations are absent, which means 
that the system is not in the OPTTR. Such a specific 
dependence of the degree of polarization of the elec- 
tromagnetic wave that passed through the sample on 
the electric field may also be a characteristic feature 
of the OPTTR. 

6. Conclusions 

To summarize, the calculations of the stationary 
and high-frequency characteristics of compensated 
GaN are carried out, which were aimed at reveal- 
ing the typical features of the streaming effect and 
the conditions needed for the effect to emerge. In 
particular, it is found that a strongly anisotropic 
distribution of electrons appears in GaN with an 
impurity concentration of 10^^ cm""^ and an elec- 
tron concentration of 10^'"' cm^'' in the range of 
applied electric fields 3 8 kV/cm and the temper- 
ature interval 30 150 K. Such a distribution mani- 
fests itself as a characteristic saturation in the de- 
pendences of the drift velocity and the total av- 
erage energy of electrons on the field. The de- 
pendence of the transverse difl'usion coefficient on 
the field decreases until it reaches the characteris- 
tic saturation. In the; framework of the low-signal 
response approximation, the spectra of the high- 
frequency electron mobility are obtained in the par- 
allel and perpendicular configurations of the station- 
ary and high-frequency fields. It is shown that, 
in the case of the parallel configuration, there ex- 
ists a transit-time resonance effect in the frequency 
range 0.5-2 THz and the field range 2-10 kV/cm, 
and the negative dynamic differential mobility can 
arise. In the perpendicular configuration, the nega- 
tive dynamic differential mobility does not arise, and 
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the spectrum of the dynamic mobility is close to the 

Drudc Lorcntzian shape. 

On the basis of aforementioned calculations, a the- 
ory is developed for the transmission of terahertz radi- 
ation through a structure with an epitaxial GaN layer. 
The relative coefBcicnt of terahertz radiation amplifi- 
cation by the structure operating in the transit-time 
resonance regime is calculated. For a single passage of 
the wave through the GaN structure, the maxinnini of 
the relative amplification coefficient is equal to 1.5%, 
which is three times as large as the losses obtaincid 
at the reflection of the same wave from metallic mir- 
rors. It is shown that the anisotropy of the dynamic 
mobility leads to a dependence of the transmission 
coefficients on the incident wave polarization. The 
polarization degree of the wave that passed through 
the structure can be controlled by changing the mag- 
nitude of applied electric field. 
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r.I. CummecbKa, B.B. Kopomeee 

EJIEKTPHHHI TA BHCOKOHACTOTHI BJIACTHBOCTI 
KOMnEHCOBAHOrO GaN B YMOBAX 
EJIEKTPOHHOrO CTPHMIHry 

P e 3 K) M e 

IIpoBefleHo fleTaJiBHHH aHajiia yinoB icnyBaHHa cTpHMlnry i 
e4)eKTa npojitoTHoro pe30HaHcy na onTHHHHx <J)OHOHax y kom- 
neHcoBaHOMy o6'eMHOMy GaN. MeTOflOM MoHTe-Kapjio npoBe- 
fleHO po3paxyHKH BHCOKOiacTOTHOi flH4'cpcimia.nbnoi pyxjiH- 
Bocxi. noKasaHO, mo npH hhbbkhx Teiniicparypax I'paTKH 30— 
77 K B ejiCKTpiiMHHx nojiHx 3—10 kB/cm b TeparepLi,OBOMy ^i- 
anaaoHi MacTOT mojkb IcnyBaTH ^HHaMlMna Bi;],'eMHa /i,n4)epeH- 
i;iajibHa pyxjmBlcTB. BnaBjieHl HOBi o3HaKH ecJieKTy cTpHMin- 
ry — aHl30Tponla flHHaMiHHoi flH4)epeHi;iajibHoi' pyxjiHBocTi 1 
oco6jiHBa noBefliHKa Koe4>iilieHTa jm^ysii y nepneHflHKyjiap- 
iiOMy nocTifiHOro ejieKxpiiMHOro nojia HanpaMKy. no6yflOBa- 
HO TBopiio npoxo;],>KeHHa TeparepLi,OBoro BHnpoMiniOBaHHa "^e- 
pe3 CTpyKTypy 3 enlTaicciftHHM mapoM GaN. OTpHMano yMOBH 
niflCHjieHHa ejieKTpoMaPHiTHHx xBHjib b fliaxiasoHi lacTOT 0,5—2 
TPn;. B ejieKTpHHHHx nojiax, 6ijibmHx, hIjk 1 kB/cm, cnocTe- 
piraeTbca nonapHsaiiiiiHa sajiejKHicTB Koe(|)l^ieHTa npoxoflJKe- 
HHa BHnpoMiHiOBaHHa Mepea CTpyKxypy. 
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